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ABSTRACT 

RX J1140.1+0307 (hereafter RX1140) is a Narrow Line Seyfert 1 (NLSl) with one of 
the lowest black hole masses known in an AGN (M < 10®Mq). We show results from 
two new XMM-Newton observations, showing soft 2-10 keV spectra, a strong excess 
at lower energies, and fast X-ray variability as is typical of this class. The soft excess 
can be equally well fit by either low temperature Comptonisation or highly smeared, 
ionised reflection models, bnt we use a covariance analysis of the fast X-ray variability 
as well as lag and coherence spectra to show that the low temperature Comptonisation 
model gives a better description of the break in variability properties between soft and 
hard X-rays. Both models also require an additional component at the softest energies, 
as expected from the accretion disc. However, this inner disc spectrum does not join 
smoothly onto the variable optical and far UV emission (which should be produced 
in the outer disc) unless the mass is underestimated by an order of magnitude. The 
variable optical and far UV emission instead suggests that L/LEdd ~ 10 through the 
outer disc, in which case advection and/or wind losses are required to explain the 
observed broadband spectral energy distribution. However, the similarity of the X-ray 
properties of RX1140 to other simple NLSl such as PG 1244+026, RE J1034+396 and 
RX J0136 means it is likely that these are also super-Eddington sources. This means 
their spectral energy distribution cannot be used to determine black hole spin despite 
appearing to be disc dominated. It also means that the accretion geometry close to 
the black hole is unlikely to be a flat disc as assumed in the new X-ray reverberation 
mapping techniques. 
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1 INTRODUCTION 

Intermediate mass black holes (IMBH), which are taken to 
refer to masses from 10“* — 10® M©, are of particular interest 
because this population would fall within the gap between 
the two principal black hole populations, i.e. super-massive 
black holes (SMBH) with M (S 10® Mq and stellar mass 
black holes with M ~ IOMq. Currently there are very few 
well studied IMBH, two examples are NGC 4395 (Filippenko 
& Ho 2003) and POX 52 (Kunth, Sargent & Bothun 1987; 
Barth et al. 2004). 

For over two decades, the most robust technique to mea¬ 
sure black hole masses in Active Galactic Nuclei (AGN) has 
been the reverberation mapping (RM) technique, which es¬ 
timates the radius of the broad line region (BLR) by mea¬ 
suring the time-lag between the line emission from the BLR 
and the intrinsic continuum emission from the accretion disc. 
Since optical RM requires well sampled spectroscopical mon¬ 


itoring of variable AGN for months or years, the total num¬ 
ber of sources with RM masses is currently less than forty 
(e.g. Kaspi et al. 2000; Peterson et al. 2004, 2005; Denney et 
al. 2006, Bentz et al. 2009b; Du et al. 2014). Instead, the RM 
measurements can be used to derive a radius-luminosity re¬ 
lation which can then give black hole mass from single epoch 
measurements using the Balmer line width (normally that of 
H/1 FWHM) and luminosity (generally monochromatic lu¬ 
minosity at 5100A, hereafter ALsioo e.g. Bentz et al. 2006, 
2009a; Denney et al. 2010). Employing this method, Greene 
& Ho (2004) identified 19 IMBH in the SDSS database. Later 
on Dong et al. (2007) and Greene & Ho (2007) extended the 
IMBH sample to include several hundred sources. 

RX J1140.1+0307 (hereafter RX1140), also referred to 
as GH 08 or SDSS J114008.71+030711.4, is amongst the 
original 19 IMBH sample of Greene & Ho (2004). It has 
a redshift of 0.081 (co-moving distance of 336 Mpc). The 
SDSS and HST images show a resolved disc component. 
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plus a bar and AGN components (Greene, Ho & Barth 
2008). The H/3 FWHM measured from the SDSS spectrum 
is 700 — 780 km (Ai et al. 2011; Jin et al. 2012a). How¬ 
ever, to get mass also requires AL 5100 , and this is made 
more difficult by the strong host galaxy contamination. De- 
convolution of the HST image gives an estimated AGN 
continuum of AL 5100 ~ 6 x lO'^'^erg s~^, which is some¬ 
what smaller than the observed variable flux of Aisioo ~ 

8 X l(Werg (Rafter et al. 2011). This gives a black hole 
mass of 8 — 10 X 10®M© (Greene, Ho & Barth 2008; Ai et al. 
2011). However, Rafter et al. (2011) applied the RM tech¬ 
nique to RX1140, and reported an upper limit of 6 light days 
for Rblr, which gives M < 5.8x IO^Mq, but this is probably 
too strict a limit given that their average sampling time of 
6.1 days. A more conservative upper limit from the reverber¬ 
ation is probably a factor 2 larger, giving M < 1.2 x 10® Mq, 
consistent with the Hf3 mass estimate. 

For this small a black hole mass, the total AGN lu¬ 
minosity is likely to be a high fraction of Eddington, with 
Zhang & Wang (2006) estimating L/LEdd ~ 1.6 from Lboi = 

9 X AL 5100 (though using the luminosity corrected for host 
galaxy contamination reduces it to ^ 0.6). So far, such high 
Eddington fraction accretion flows, with L/LEdd 1, are 
found in the subset of broad line AGN known as Narrow- 
Line Seyfert 1 (e.g. Osterbrock & Pogge 1985; Boroson & 
Green 1992; Leighly 1999). Their narrow Balmer lines in¬ 
dicate relatively small black hole masses ( 1 O®“^M 0 : Boro¬ 
son 2002 ), so that their bolometric luminosities are close to 
Eddington. For such high Eddington fractions it is clearly 
a moot point as to whether the H/3 scaling relations hold. 
Marconi et al. (2009) suggest a correction to the mass scal¬ 
ings such that the BLR clouds trace the effective gravity, 
but the lack of observed wind features in the low ionisation 
BLR lines form a strong argument against such effects being 
important. This could indicate that the clouds are optically 
thick to electron scattering, with columns of > 10 ^"^ cm“^, 
or that the radiative acceleration only affects the front face 
of the cloud (Baskin, Laor & Stern 2014). Thus it seems 
probable that the H/3 line widths reflect the black hole mass 
even at Eddington and potentially beyond. 

One of the ubiquitous features of high Eddington frac¬ 
tion AGN is their prominent soft X-ray excess superposed 
on a steep X-ray power law which is very variable. These 
properties are seen in the X-ray spectrum of RX1140 (Mini- 
utti et al. 2009; Ai et al. 2011). While the variability scales 
with black hole mass, making it unsurprising that the X-ray 
variability of RX1140 is amongst the strongest seen AGN 
(e.g. Ponti et al. 2012; Ludlam et al. 2015), the origin of the 
soft excess is more puzzling. Models which reproduce the 
shape include optically thick, low temperature Comptonisa- 
tion (e.g. Laor et al. 1997; Magdziarz et al. 1998; Gierlihski 
& Done 2004; Jin et al. 2013), ionised, highly smeared reflec¬ 
tion (e.g. Fabian & Miniutti 2005; Zoghbi et al. 2010; Fabian 
et al. 2013; Fabian, Kara & Parker 2014; Uttley et al. 2014), 
and smeared absorption (Gierlihski & Done 2004; 2006). 
The smeared absorption model is disfavoured on theoreti¬ 
cal (Schurch & Done 2007) as well as observational grounds 
(Miniutti et al. 2009, Ai et al. 2011), but both Gomptoni- 
sation and smeared reflection provide a good fit to the soft 
excess in this object (Ai et al. 2011). 

However, for such low mass, high mass accretion rate 
AGN, the disc itself should also contribute a significant 


fraction of the emission in the soft X-ray bandpass. It is 
well known that the accretion flow spectrum is not well 
described by a disc in standard broad line Seyfert Is and 
Quasars (e.g. Elvis et al. 1994), but the higher Eddington 
ratio NLSl have spectra which do appear to be more disc- 
dominated (Jin et al. 2012a,c; hereafter J12a; c; Done et al. 
2012, hereafter D12). This is most easily seen in the context 
of the full Spectral Energy Distribution (SED) as the opti- 
cal/UV disc normalisation is determined by the product of 
mass and mass accretion rate, with L^ oc Thus 

this directly measures mass accretion rate if the black hole 
mass is known, and hence the total luminosity L = rjMc^, 
where the efficiency rj depends on black hole spin (e.g. Davis 
& Laor 2011; Done et al. 2013). This also predicts where the 
disc disc emission peaks, and show that the disc spectrum 
should extend into the soft X-ray bandpass for the NLSl 
even at low spin (Done et al. 2013). 

RX1140 is the lowest mass NLSl included in the sam¬ 
ple of 51 unobscured Type 1 AGN of J12a,c. These papers 
carried out a detailed study of the broadband SED from 
optical to hard X-rays, fitting them with the new accretion 
models which allow some of the accretion energy to be dis¬ 
sipated in a Gomptonised soft X-ray excess and power law 
tail as well as in the thin disc (D12). These assume energy 
conservation i.e. that there is no strong energy loss via ad- 
vection and/or winds as might be expected in highly super- 
Eddington sources. J12a, c used the black hole mass was a fit 
parameter within the uncertainties derived from using the 
intermediate and broad H/3 line component widths as lower 
and upper limits, respectively. This gave a best fit value 
of 2.9xlO®M0 for this source (J12a), which increased to 
6x 1O®M0 when the updated SED models including a colour 
temperature correction to the disc emission were used (D12, 
J12c). Clearly the mass of RX1140 derived from broadband 
SED fitting is much higher than the virial and RM mass 
estimates, and this requires further consideration. 

The X-ray properties of RX1140 are very similar to 
a small group of the most extreme NLSls, such as PG 
1244+026, RE J1034+396 and RX J0136.9-3510, which all 
exhibit strong X-ray variability, prominent and featureless 
soft excesses and are all reported to be accreting around the 
Eddington limit (Middleton et al. 2009; Jin et al. 2009; Jin 
et al. 2013). In order to investigate the physical interpreta¬ 
tion for the soft excess seen in RX1140, a detailed timing 
analysis must be performed in conjunction with the spectral 
analysis (e.g. Jin et al. 2013). In achieve this we applied for 
a long XMM-Newton observation for this source. This was 
successful, but was split into two observations separated by 
two weeks. In this paper we present the results from a com¬ 
bined analysis of these two newXMM-Newton observations, 
plus the previous one extracted from the archive. 

This paper is organised as follows: Section [2] describes 
the data reduction procedures used for the XMM-Newton 
observations. Section [ 3 ] presents variability properties of 
RX1140. Section [ 4 ] presents a combined spectral and vari¬ 
ability study of the target. We compare two physical mod¬ 
els for explaining the properties of the 0.3-10 keV spectra 
and variability. In Section^ we perform broadband SED fit¬ 
ting by including ultraviolet (UV), optical and near-infrared 
data. In section[^we discuss issues such as various estimates 
of the black hole mass, UV/optical luminosity and the mech¬ 
anism responsible for the soft X-ray excess. Our summary 
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and conclusions are given in section [T] When converting 
flux to luminosity, we adopt a flat universe model with the 
Hubble constant Ho = 72 km s~^ Mpc~^, Q.m ~ 0.27 and 
= 0.73. 


2 DATA REDUCTION 

RX1140 was observed by the XMM-Newton satellite on 3rd, 
Dec., 2005 (Principal Investigator: Dr. Giovanni Miniutti) 
(hereafter: Obs-1). We proposed a longer observation with 
the aim to study its X-ray variability. Our observation was 
divided into two parts due to satellite scheduling. These were 
carried out on 18th, Dec., 2013 (hereafter: Obs-2) and 1st, 
Jan., 2014 (hereafter: Obs-3). During the three observations, 
all EPIC cameras were in the full window mode. We used 
SAS vl3.5.0 and the latest calibration files, and followed the 
standard procedures to reduce the data. We chose the source 
extraction region to be a circular region of radius 45”for each 
EPIC camera. The background was selected from a nearby 
circular region with the same radius as for the source. For 
Obs-1 the net source count rates are 0.62 ct s~^, 0.13 ct s~^ 
and 0.13 ct s~^ for the PN, MOSl and MOS2, respectively. 
Obs-2 has a mean count rate that is half that of Obs-1. Obs- 

3 has a similar mean count rate as Obs-1. All of these count 
rates are well below the threshold count rates capable of 
causing a photon pile-up effect in the full window mode of 
each EPIC camera. 

We selected data with PATTERN 12 for MOSl and 
MOS2, and PATTERN ^ 4 for the PN. Light curves were ex¬ 
tracted from both the source and background regions. There 
were high background flares in all three observations. There¬ 
fore we visually checked the background light curves to iden¬ 
tify high background sections and removed them. Then the 
background was subtracted from the source light curve us¬ 
ing LCMATH in FTOOLS. Spectra were extracted for the source 
and background regions, separately. Response matrices were 
produced using RMFGEN and ARFGEN. Areas for the source 
and background regions were calculated using BACKSCALE. 
We also used the SAS task RGSPRDC to extract the 1st order 
spectra from RGSl and RGS2. Response matrices for the 
RGS were produced with RGSRMFGEN. All spectra were re¬ 
binned by GRPPHA with a minimum of 25 counts per bin, so 
that the fitting is appropriate.But the signal-to-noise of 
RGS spectra are too low to provide further information, so 
we do not present them in this paper. All spectral fittings 
were performed in xspec v12.8.2 (Arnaud 1996). 

There are also simultaneous Optical/UV observations 
obtained using the XMM-Newton Optical Monitor (OM). 
OBs-1 has UVWl and UVM2 filter data, while Obs-2 and 
3 have U, B and UVWl. We searched the OM source list 
file to obtain the count rate of RX1140, in each available 
Alter, and inserted these values into the standard OM data 
file template om_filter_default.pi. The data file was then com¬ 
bined with the ‘canned’ response filefl to be ready for xspec 
fitting. 


^ http://heasarc.gsfc.nasa.gov/FTP/xmm/data/responses/om 
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Figure 1. The 100 s binned, background subtracted, 0.3-10 keV 
light curves for Obs-1, Obs-2 and Obs-3. Every cyan region is a 
11 ks segment that is free from high background contaminations. 


3 VARIABILITY STUDY 

The background subtracted PN light curves over the total 
energy band (0.3-10 keV) are shown in Figure [T] Strong 
variability is observed over all timescales. This is also con¬ 
firmed by their Power Spectral Density (PSD) in Figure [S] 
The flat PSD indicate that intrinsic variability (i.e. Poisson 
noise subtracted) exists throughout the 0.2-10 ks timescale. 
None of the observations show any sign of a high frequency 
break. This is clearly different to the case of PC 1244-1-026, 
where the 0.3-1 keV power spectrum drops substantially 
above ~ 10“® Hz (Jin et al. 2013), thus strongly suggests 
that RX1140 is substantially lower mass than PC 1244-1-026, 
and that the mass is < 10 ®Mq (Ponti et al. 2012; Kelly et 
al. 2013; Ludlam et al. 2015). 

The time-averaged X-ray spectra for Obs-1 (black) and 
Obs-2 (red) are shown in the upper panel of Figure |3l The 
mean spectrum of Obs-3 is almost identical to Obs-1, so it 
is not shown in the figure. The Obs-1 and Obs-2 spectra are 
quite similar, but a ratio of the best fit model to Obs-1 (black 
solid model), rescaled to the observed 2-10 keV spectrum of 
Obs-2 (red solid model) under-predicts the soft emission. 
This longer term (weeks to years) variability can be com¬ 
pared to the rapid variability. The shapes of high [pnjrate ^ 
0.7 ct s~^: orange) and low {pnjrate ^ 0.35 ct s~^: blue) 
count rate spectra from Obs-1 are identical within errors to 
the mean spectrum. This can be seen in the lower panel of 


























4 C. Jin, C. Done, M. Ward 



Figure 2. Power Spectra Density for the 0.3-10 keV lightcurves of the three observations obtained using FTOOL powspec after subtracting 
the Poisson counting noise. For Obs-3, only the two 11 ks segments in Figure[T]were used to avoid the influence of high background. 


Figure O which gives the ratio of each dataset to the mean 
spectral model renormalised to the 2-10 keV count rate. 


3.1 RMS Spectra 

Fractional variability observed in different energy bins, i.e. 
the RMS spectra (e.g. Edelson et al. 2002; Markowitz, Edel- 
son & Vaughan 2003; Vaughan et al. 2003), can be very use¬ 
ful in decomposing the different spectral components. We 
use the segment algorithm detailed in Wilkinson & Uttley 
(2009) to calculate the RMS. For every energy bin, the light 
curve is divided evenly into several segments. Then the ex¬ 
cess variance is calculated for every segment and averaged 
over all segments. The square-root of the excess variance 
is RMS. The error-bars are calculated using Equation B2 in 
Vaughan et al. (2003). The binning time and segment length 
determine the timescale for which the RMS is calculated. 
The high background intervals interrupt the continuity of 
the light curve. In Obs-3, the longest continuously sampled 
segment is only 11 ks, while it is longer in Obs-1 and Obs-2. 
To mitigate bias from high background contamination, we 
choose a longest segment length of 11 ks, and so are able 
to identify a total of eight segments with little high back¬ 
ground, i.e. three segments in Obs-1, three in Obs-2 and two 
in Obs-3 (see the cyan regions in Figure [T]). 

There are insufficient data to produce frequency- 
resolved RMS spectra. Instead, we calculate only RMS of the 
highest frequency (HF) variability using timescales of 0.2-2 
ks from the eight llks segments. The results are shown in 
Figure m Interestingly, the HF RMS spectra reveal a rising 
shape from soft to hard X-ray ranges in Obs-1 and Obs-3, 
which is similar to that observed in some other high mass ac¬ 
cretion rate NLSls, such as PG 1244-1-026 (Jin et al. 2013), 
RE J1034-b396 (Middleton et al. 2009) and RX J0136.9-3510 
(Jin et al. 2009). This difference in variance between soft and 
hard energy bands can be easily explained if the soft X-ray 
excess varies less than the power law tail, favouring mod¬ 
els where this is a true additional component e.g. from soft 
Comptonisation (e.g. Gierliiiski & Done 2004, Middleton et 
al. 2009, Jin et al. 2013, Matt et al. 2014) rather than due 
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Figure 3. Time-averaged spectra of Obs-1 (black) and Obs-2 
(red) together with the highest (orange) and lowest (blue) spectra 
seen within Obs 1. The best-fit model of Obs-1 is the same as in 
Figure [5l (ai (black model), then renormalised to match the 2-10 
keV of each spectrum. The lower panel shows the ratio of each 
dataset to the model, showing significant spectral variability only 
in Obs-2, where the soft X-rays are weaker. 


to ionised, smeared reflection, as this tends to contribute 
equally in soft and hard bands (Gierlitiski & Done 2006). 

Another key feature seen in Figure 0] is that the 0.3-1 
keV HF RMS in Obs-2 is rather different to that of Obs-1 
and Obs-3, with much less difference between the fast vari¬ 
ability as a function of energy (more HF power below 1 
keV, and marginally less at the highest energies). Whatever 
happened to depress the soft spectrum relative to the mean 
(Figure O) may also increased its fast variability (see Gard¬ 
ner & Done 2015 for a possible explanation of this in terms 
of cloud occultions in the inner disc). 
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Figure 4. The HF (0.2-2 ks) RMS spectra for the three observations. Obs-1 and Obs-3 have similar fast spectral variability as well as 
similar mean spectra, while Obs-2 has more fast variability at low energies than both Obs-1 and Obs-3. 


3.2 Covariance spectra 

The RMS variability spectrum shows all variability at each 
energy, but where there may be multiple components. A 
more powerful technique is the covariance spectrum (devel¬ 
oped by Wilkinson & Uttley 2009) as it determines the spec¬ 
trum of the variability correlated with a given energy band 
lightcurve used as a reference. For the NLSl PG 1244-1-026, 
Jin et al. (2013) showed that using the fastest variability 
in the hard X-ray (2-10 keV) lightcurve as a reference band 
gave strong evidence for a Comptonisation origin for the soft 
X-ray excess, while using the fastest variability in the 0.3- 
1 keV band showed evidence for a contribution at the lowest 
energies from the disc. 

Here we apply the same technique to our datasets to 
produce the HF (0.2-2 ks) covariance spectra using hard (2- 
10 keV) and soft (0.3-1 keV) X-rays as the reference band. 
We use only PN data so as not to have differences in response 
from combining PN and MOS. But the data quality means 
that the signal-to-noise of the covariance spectra for RX1140 
is much lower than for PG 1244-1-026. Then we study both 
the covariance spectra and the time-averaged spectrum (see 
Figure [5]) in order to try to determine the origin of the soft 
X-ray excess, i.e. whether it can be better fit by an addi¬ 
tional Comptonisation model or by relativistically smeared 
reflection, and whether there is also a contribution from the 
accretion disc itself. 


4 SPECTRAL ANALYSIS 

The X-ray spectrum derived from Obs-1 has been analysed 
previously. Miniutti et al. (2009) found that an X-ray reflec¬ 
tion plus a thermal disc model, gives a better fit than the 
smeared absorption model. Ai et al. (2011) found that re¬ 
flection and Comptonisation models gave comparably good 
fits, and that both these were better than smeared absorp¬ 
tion and p-free disk models. The shape of the HF RMS and 
mean spectra in Figure |3] and Figure U also indicates that 
RX1140 may be similar to other extreme NLSls such as 
PG 1244-f026, RE J1034-f396 and RX J0136.9-3510, whose 


soft X-ray excess can be well interpreted by Comptonisa¬ 
tion model. Therefore, we try both Comptonisation and re¬ 
flection models for Obs-1 and Obs-2 spectra. The best-fit 
models are shown in Figure [5] 


4.1 Comptonisation Model 

For the Comptonisation scenario, we use a low tempera¬ 
ture, optically thick Comptonisation model (comptt model, 
Titarchuk 1994; Hua & Titarchuk, 1995; Titarchuk & 
Lyubarsky 1995) to fit the soft X-ray excess; a high temper¬ 
ature Comptonisation (nthcomp model in xspec, Zdziarski, 
Johnson & Magdziarz 1996; Zycki, Done & Smith 1999) to fit 
the power law tail, together with its neutral, smeared reflec¬ 
tion (kdblur*pexmon, Laor 1991; recoded as a convolution 
model, Nandra et al. 2007). We also include an accretion 
disc component (diskbb model, Mitsuda et al. 1984; Mak- 
ishima et al. 1986). Using the equation in Peterson (1997) 
and assuming = 6rg, Log{M) = 5.77 and L/LEdd = 2.69 
(Miniutti et al. 2009), we calculate the temperature at the 
inner radius (Tin) to be 110 eV. Including a colour tem¬ 
perature correction would increase this, while allowing the 
some of the inner disc power to be dissipated in powering the 
soft and hard X-ray Compton components would decrease 
it, so we freeze the parameter Tin at this temperature. It 
also serves as the temperature of the seed photons for both 
soft and hard X-ray Comptonisation components. We fix 
the Galactic column at Nh = 1.91x10^° cm~^ (Kalberla 
et al. 2005) but allow additional free absorption by the host 
galaxy gas column (zwabs, Morrison & McCammon 1983). 
The complete xspec model parameters are given in Table [T] 
The HF variability is likely to originate mainly from 
the inner region of the flow, where the hard X-ray emission 
is produced. Hence we expect that the spectrum correlated 
with the HF 2-10 keV variability should have the same shape 
as the high temperature Gomptonisation component which 
dominates the 2-10 keV spectrum in this model. However, 
reflection (and associated thermal reprocessing which could 
contribute to the soft X-ray excess: Gardner & Done 2015a) 
should also vary in a correlated manner as the expected 
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Figure 5. X-ray spectral fitting of Obs-1 (Panel Compt-Obsl, Refi-Obsl) and Obs-2 (Panel Compt-Obs2, Refi-Obs2). The upper row 
shows the results of Comptonisation fitting to the mean (red) and HF covariance (blue) spectra. In each panel, red solid line is the best-fit 
model, red dashed line is the hard X-ray Comptonisation component, red dotted line is the soft X-ray Comptonisation component, red 
dash-dot line is the hot accretion disc component. The blue lines are spectral fits to the HF covariance spectra with only the soft and 
hard Comptonisation components. The lower row panels present the results of reflection model fitting. The red dashed line is the intrinsic 
hard Comptonisation component, the red dotted line is the reflected component by accretion disc, the red dash-dot line is the accretion 
disc component. The blue lines are spectral fits to the covariance spectra with only the hard Comptonisation and reflection components. 


reverberation timescale (Fabian et al. 2009) is short for a 
low mass black hole. Only the disc is not expected to have 
any variability on these timescales. 

To perform spectral fitting, we first applied the spec¬ 
tral model to the entire 0.3-10 keV spectrum, but we found 
some parameters were not well constrained because of pa¬ 
rameter degeneracy. Therefore we attempt to constrain the 
hard X-ray Comptonisation and neutral reflection compo¬ 
nents by applying them only to the 2-10 keV spectra of 
Obs-1 and Obs-2. The model fits the 2-10 keV spectra 
well (xJ = 84/82). The hard X-ray photon index (T) is 
found to be 2.26±0.12, the inner radius parameter (Rin.) 
of kdblur is 20lfg'’ with reflection solid angle n/27r fixed 


at unity. Removing the neutral reflection component leads 
to r = 2.13 ± 0.12, but the fitting statistics is only slightly 
worse (xS = 96/84). The low signal-to-noise in the hard X- 
ray band cannot constrain this reflection component well. 
Then we re-fit the model to the entire 0.3-10 keV band by 
fixing r at 2.26 and Rin at 20, and all spectral components 
are better constrained then before. The best fit parameters 
are listed in Table [T] 

In Obs-1, the HF covariance spectrum can be well fit¬ 
ted by the hard Comptonisation with F = 2.26, without 
significant contribution from a variable soft Comptonisation 
component. On the contrary, the HF covariance spectrum in 
Obs-2 is softer so requires more contribution from the soft 
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Table 1. The XSPEC model and the best-fit parameters in Figure [5] The upper and lower limits are the 90% confidence range. 


Model Name 

Model Expression in XSPEC vl2.8.2 


Comptonisation 

C0NSTANT*WABS*ZWABS*( DISKBB+ 

NTHCOMP + COMPTT + 

KDBLUR*PEXM0N ) 

Spectrum 

Component 

Parameter 

Value 

Value 




Obs-1 

Obs-2 

Meein Spec 

ZWABS 

Nff (10^^ cm~'^) 

rv +0.82 
^ -0 

0 027 

-0.015 


DISKBB 

Tin (keV) 

0.11 fixed 

0.11 fixed 


DISKBB 

norm 

135 

135 fixed 


NTHCOMP 

r 

2.26 fixed 

2.26 fixed 


NTHCOMP 

kTseed (keV) 

tied to Tin 

tied to Tin 


NTHCOMP 

kTe (keV) 

100 fixed 

100 fixed 


NTHCOMP 

norm 

(1.26 lo;i°)xlO-4 

(8.11 7S:8o)xio"" 


COMPTT 

kTseed (keV) 

tied to Tin 

tied to Tin 


COMPTT 

kT (keV) 

0.38 

0 21 +1-^^ 

-0.07 


COMPTT 

r 

10.9 

19.2 71?:^ 


COMPTT 

norm 

(2.68 71;°®) X10-3 

(1.81 71:?^)xio-3 


KDBLUR 

R-in (Rp) 

20 fixed 

20 fixed 


PEXMON 


-1.0 pegged 

-1.0 pegged 

Cov. Spec 

NTHCOMP 

norm 

(5.86 

(2.21 7767)x10-5 


COMPTT 

norm 

(1.45 71;2)xl0-4 

(9.74 79^7^(xlO-'l 



xl 

298/314 = 0.95 

176/224 = 0.78 

Reflection 

C0NSTANT*WABS*ZWABS*( DISKBB 

+ NTHCOMP + KDBLUR*RFXC0NV*NTHC0MP ) 

Spectrum 

Component 

Parameter 

Value 

Value 




Obs-1 

Obs-2 

Mean Spec 

ZWABS 

Nif (10^^ cm~'^) 

+0.02 
^ -0 

u.uuo _o.oo5 


DISKBB 

Tin (keV) 

0.11 fixed 

0.11 fixed 


DISKBB 

norm 

211 

-136 

83 786’’ X 10-4 


NTHCOMP 

r 

2.42 

o r)Q +0.10 
-0.12 


NTHCOMP 

norm 

5.72 7°;®° X 10-5 

3.60 71?? X 10-5 


KDBLUR 

Index 

5 fixed 

5 fixed 


KDBLUR 

R-in (Rg) 

4.98 ttil 

3-05 hlf. 


RFXCONV 

R^lre/Z 

-2.19 7L^" 

-2.29 7L®’ 


RFXCONV 

^^abund 

1 39 +2.50 

L.ov _o.7i 

0.81 t°o:Ti 


RFXCONV 

log i 

o njY +0.37 
-0.22 

2.86 7?:?? 

Cov. Spec 

NTHCOMP 

norm 

4.90 7^5? X 10-5 

1.08 71?^ X 10-5 


RFXCONV 

R®lre/Z 

-0.28 7^ 

-2.60 71,"° 



xl 

295/312 = 0.95 

170/222 = 0.77 


Comptonisation (Figure [5^ and Figure [Dd). Neither covari¬ 
ance spectrum requires any disc component. 

We can also test for the presence of an accretion disc 
component in the soft X-ray bandpass by removing the disc 
component from the model and then re-fitting the Obs-1 
spectra. The electron temperature (kT) of the soft Comp¬ 
tonisation component increases from 0.38 keV to 0.46 keV, 
the optical depth (r) decreases from 10.9 to 7.1, the normal¬ 
isation increases by 10%. However, xS increases from 0.94 
to 0.98 (i.e. Ax^ = 10 for the removal of 2 free parame¬ 
ters, which is only marginally significant). The significance 
increases if we consider only the 0.3-0.5 keV band, where xS 
increases from 1.22 to 1.33 after removing the disc compo¬ 
nent. Thus both the time average spectra and HF variabil¬ 
ity support the decomposition of the soft X-ray spectrum 
into at least two components, which we interpret here as a 
contribution from the accretion disc, and an additional low 
temperature Comptonisation. A high temperature Compton 
spectrum is required to reproduce the 2-10 keV spectrum. 


4.2 Reflection Model 

Ionised, strongly smeared disc reflection plus a thermal 
accretion disc component was previously reported to be 
a good model for the time-averaged Obs-1 spectrum of 
RX1140 (Miniutti et al. 2009, Ai et al. 2011). Thus we use 
diskbb+nthcomp+kdblur*nthcomp as our spectral model, 
with Tin of the disc component fixed at 110 eV as before. 
Figure [SJ; and Figure [SJi shows this model applied to Obs-1 
and Obs-2 respectively. Again, the low statistics mean that 
some parameters cannot be well constrained. We choose to 
fix the illumination index of the smeared reflection at 5, and 
fit for the inner disc radius. This is always small in this 
model, at ~ 3 — 5Rg , so any reverberation timescale should 
be short. Hence both the high temperature Comptonisation 
and its reflection should vary together. However, the rela¬ 
tively flat HF covariance spectrum in Obs-1 favours a less 
reflection component contribution than in the time averaged 
spectrum, while the softer HF covariance spectrum in Obs-2 
favours a larger reflection contribution. But again the poor 
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Figure 6. Frequency and energy resolved time-lag and coherence for RX1140. The points are derived from five segments of 11 ks length 
each, which comprise three segments in Obs-1 and two segments in Obs-3 (see Section Ell. 


signal-to-noise means that this is not very significant (Ta¬ 
ble [T]). Similar as in the Comptonisation scenario, accretion 
disc component is not detected in the HF covariance spectra. 


4.3 Lag and coherence spectra 

Both Comptonisation and reflection models for the soft X- 
ray excess can give comparably good fits to both the time 
averaged and HF covariance spectra. Both models also re¬ 
quire that there is an additional soft component, which we 
interpret as the direct disc emission, at the softest energies. 

We can also use spectral coherence and time lags to 
constrain the X-ray spectral mechanisms (e.g. IvT^Hardy et 
al. 2004; Fabian et al. 2009). However, the low count rate 
of RX1140 and the relatively short continuously sampled 
light curve segments (<11 ks, see Section [Sj mean that the 
statistics are limited. We calculate these as in Nowak et al. 
(1999), using only Obs-1 and 3 (since Obs-2 has a signif¬ 
icantly different spectrum and variability). We choose low 
and high energy band with 0.3-0.7 keV and 1.8-2.5 keV so 
as to try to best isolate different spectral components whilst 
still having enough counts. 

Figure [6^ shows that there is a soft reverberation lag 
of ~ 200 s seen between the 0.3-0.7 keV and 1.8-2.5 keV 
light curves at a frequency of ~ 2 x 10“® Hz. But the zero 
spectral coherence and the dominance of Poisson noise at 
this frequency means this lag can be artificial. However, 


there is a clear positive lag at the lowest frequencies (be¬ 
low 2 X 10“"* Hz), where the soft leads the harder band. 
This lag is supported by the good spectral coherence be¬ 
tween these bands up to 6 x 10“'* Hz (Figure [6j: and Figure 
Ell). Note that Poisson noise will not be an issue on such 
long timescales. Figure [Sli) shows the energy spectrum of 
this long timescale soft lead, showing that there is a clear 
break between the time lags at low and high energy, although 
we do not have enough statistics at hard X-ray. This is not 
easily compatible with the reflection model shown in Fig¬ 
ure Eb), where a single component dominates the spectrum 
from 0.4-10 keV. 

Thus while the signal-to-noise is poor, these data pro¬ 
vide marginal evidences in terms of spectral coherence and 
time-lag to support a model where the majority of the soft 
X-ray excess does not contribute to the 2-10 keV bandpass, 
favouring the low temperature Comptonisation origin over 
the reflection dominated model. 

5 BROAD BAND SPECTRAL ENERGY 
DISTRIBUTION 

5.1 Multi-wavelength Data 

As outlined in the introduction, RX1140 was reported as an 
IMBH, whose mass is estimated to be 5 —10 x 10®M© accret¬ 
ing at a high fraction of the Eddington limit. We construct 
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a multi-wavelength SED based on non-simultaneous obser¬ 
vations (Figure [T]). We use the XMM-Newton PN spectrum 
from Obs-1 taken in 2005, together with archival data from 
the ROSAT PSPCB (1992-06-16) which we reduced using 
XSELECT v2.4b. The spectra were regrouped fixing a mini¬ 
mum of 25 counts per bin, and are in good agreement with 
the time-averaged PN XMM-Newton spectrum. 

In the optical and UV band, there are simultaneous pho¬ 
tometric points from XMM-Newton OM for Obs-1 (UVM2 
and UVWl filters) and Obs-2 and Obs-3 (UVWl, U and B 
filters). The UVWl filter fluxes of Obs-1 and Obs-3 are very 
similar, so are their X-ray fluxes, so we use all OM data from 
these two observations (except for UVWl since its waveband 
is already covered by the combination of UVM2, U and B 
filter bands). We do not use Obs-2 as it is clearly different in 
optical/UV as well as X-ray, with only 90% flux in UVWl, 
84% in U and 99% in B compared to Obs-3. These are less 
suppressed than the X-ray flux (Obs2/Obs3 = 46% in 0.3- 
2 keV, 63% in 2-10 keV), indicating that the amplitude of 
UV/optical variability is less than that of the X-ray over 
variability timescales of weeks and years. 

Other data points in Figure [7] include GALEX pho¬ 
tometry (obs-date: 2007-03-22), SDSS UGRiz photometry 
and fibre spectrum (obs-date: 2001-03-26), UKIDSS Pet¬ 
rosian magnitude of JYHK photometry (obs-date: 2008-04- 
21). The difference in the normalisation, especially that for 
the UKIDSS and SDSS fibre fluxes, is likely due to different 
aperture sizes. 

5.2 Broadband SED Fitting 

We use optxagnf in xspec v12.8.2, to fit the SED from the 
optical to hard X-ray energies (D12). This includes both 
intrinsic disc emission, together with low and high temper¬ 
ature Comptonisation, as favoured by the spectral fits in 
Sectional An additional neutral reflection (kdblur*pexmon 
model, Laor 1991, recoded as a convolution model by Nandra 
et al. 2007) is added to fit the marginally rising shape above 
4 keV. We fix the inclination of the reflector at 30°, and fix 
the spectral index of the hard Comptonisation component 
at 2.26 (see Section[4|. We correct for extinction by a Galac¬ 
tic column density of Nh = 1.91x10^° cm~^ (Kalberla et 
al. 2005). We also allow for intrinsic extinction in the host 
galaxy. This gives a best fit around Nh = 2x10^° cm~^ 
at z = 0.0811 (zwabs model, Morrison & McCammon 
1983). Both Galactic and intrinsic reddening are assumed 
to follow the standard dust to gas conversion formula of 
E{B — V) = 1.7x 10“^^ Vir (Bessell 1991) and modelled by 
redden and zredden model (Cardelli et al. 1989) in XSPEC. 

The far UV GALEX points fit nicely on the expected 
disc spectrum from optxagnf, but the optical and near 
UV points display a curvature similar to that of host galaxy 
emission. The SDSS images show that the host galaxy of 
RX1140 is of an Sc type galaxy. Therefore we select an Sc 
galaxy template from the SWIRE library (Polletta et al. 
2007) and incorporate this into our XSPEC model. A constant 
normalisation offset is allowed between the UKIDSS points 
and other data. 

Figure [7| shows the extinction/reddening/aperture cor¬ 
rected data (coloured points, compared to the observed data 
points in grey) and the best-fit SED model. The formal fit¬ 
ting statistics are poor (reduced ~ 982/330), which is 



Figure 7. Broadband SED of RX1140. Black data in the X-ray 
band is the Obs-1 mean spectrum, purple data in the soft X- 
ray band is from ROSAT. The UV, optical, infrared photometric 
points comprises GALEX (blue circles), OM of XMM-Newton 
(purple diamonds), SDSS UGRIZ (green squares) and fibre (green 
spectrum), UKIDSS JYHK (cyan stars). The best fit optxagnf 
model assuming black hole mass as a free parameter is shown as 
a red solid line is the best-fit model, orange dotted lines show 
different SED components. This gives M ~ IO^Mq, much higher 
than expected from H/3 line width and reverberation mapping. 


mainly due to the small error bars of the photometric points 
in the optical/UV band. 

However, the best-fit black hole mass is 9.6 x 10®M© 
(with correspondingly fairly low L/LEdd = 0.17). This is al¬ 
most an order of magnitude higher than the expected mass 
from optical spectroscopy and reverberation (Greene, Ho & 
Barth 2008; Rafter et al. 2011). The Eddington ratio is cor¬ 
respondingly low, L/LEdd = 0.17, which is in conflict with 
the classification of this object as a NLSl, as supported by 
the steep X-ray spectrum, and by the optxagnf coronal ra¬ 
dius of ~ 15i7g, consistent with the typical value for NLSl 
(J12c, D12). 

We test how robustly the SED modelling constrains 
mass by fitting pure disc models down to 15779 (dashed 
lines) to the UV points for masses of 1.5 x 10® (blue: 
L/LEdd = 400), 1.0 X 10® Mq (orange: L/LEdd = 10) and 
1.0 X 10^ Mq (red: L/LEdd = 0.17 as before. The red solid 
line shows the previous best fit. These disc component all 
give optical fluxes at SIOOA which match to the observed 
variable flux (Rafter et al. 2011). However, all the reason¬ 
able masses over-predict the observed X-ray emission by fac¬ 
tors of 10-100. Allowing the disc to extend down to &Rg, or 
1.2Rg for a maximally spinning black hole, would make this 
discrepancy worse (e.g. Done et al. 2013). 

In the right panel of Figure [S] we use the same range of 
masses but aim instead to primarily fit the X-ray spectrum, 
and extrapolate the model to lower energies. The model with 
a 10^ Mq (red curve) can fit the far UV (and variable opti¬ 
cal) data, whereas the model for the 10® and 1.5 x 10® Mq 
black hole masses under-predict the disc continuum by a fac¬ 
tor of 10-50 as the peak of the disc spectrum shifts into the 
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Figure 8. Broadband SED fitting similar as in Figure[7l but with different black hole masses: M = 9.6 x IO^Mq (red), 1.0 x IO^Mq 
(blue), 1.5 X 10^ Mq (green). In Panel-a, the models primarily fit the UV points and have L/L^dd = 0.17,10 and 400, respectively for disc 
emission outside 15 Rg (see Section E21i. In Panel-b, the models tend to primarily fit the X-ray spectra and have LjLEdd = 0.17,0.56 
and 2, respectively. Solid lines are the total models, while dashed lines are the disc components. 


soft X-ray band, and its contribution to the UV and optical 
decreases. 

6 DISCUSSION 
6.1 Black Hole Mass 

Section [T] summarised previous mass estimates for RX1140 
from RM (Rafter et al. 2011) and H/3. All previous black 
hole masses based on optical Spectra report RX1140 to be 
an IMBH with M < 10® Mq. 

The X-ray variability power spectra also strongly sup¬ 
port such a low mass, with no observed high frequency 
break on timescales longer than 500 s (see Figure. [2|. This 
means that the total excess variance, arms, is high. We 
include our new data on the M — a^ms relation (Mini- 
utti et al. 2009; Zhou et al. 2010; Ponti et al. 2012). Fig¬ 
ure [9] shows the reverberation mapped subsample of ob¬ 
jects from Ponti et al. (2012), where a^ms values were cal¬ 
culated from 40 ks light curves in the 2-10 keV bandpass. 
The best fit linear regression line is also shown. For RX1140, 
arms = 0.11±0.04 in Obs-1 and 0.12±0.08 in Obs-2 (Obs-3 
does not have sufficiently good data). These values corre¬ 
spond to M — 1.25xlO®M0. A Mq mass lies outside 
the 2a region, but any mass within 3x10® — 3.0 x10®Mq is 
consistent with the relation (see also Ludlam et al. 2015). 

This highlights the inconsistency between the black hole 
mass derived from the SED models of 10^Mq, and all 
other mass determination methods which give 10® Mq. For 
a mass of 10® M©, the SED models can either fit the opti- 
cal/UV flux from the AGN, but then over-predict the X- 
rays, or they can fit the X-rays, but under-predict the opti- 
cal/UV (Figure ISJr). Given that the fraction of optical/UV 
flux from the AGN is constrained by the observed variabil¬ 
ity (Rafter et al. 2011), then this requires that the 10®M© 
disc fit in Figure [8^) correctly describes the outer disc, with 


a derived L/LEdd = 10. Clearly this is substantially super- 
Eddington, so there can be substantial energy loss via op¬ 
tically thick advection and/or winds from the inner disc. 
The fitted SED models, which assume that the disc is sub- 
Eddington so that neither of these processes is important 
(orange and blue curves in Figure fig-sed2a), then over¬ 
estimate the total luminosity. 

The extent of the energy losses due to winds/advection 
can be estimated from the best fit SED model. The pa¬ 
rameters of M = 10^Mq,L/LE dd = 0.17 give Ltoi ~ 
2 X 10"^"^ ergs cm~^ s“^, which is around the Eddington limit 
for the most probable mass of 10® M©. This is a supporting 
evidence for the proposal of Wang et al. (2014) that super- 
Eddington black holes have total luminosity which saturates 
around the Eddington limit. 

However, RX1140 also has another unusual feature in 
that it exhibits strong optical variability, as seen in the RM 
campaign of Rafter et al. (2011). The majority of NLSl show 
less optical variability than the broader line Seyferts (e.g. Ai 
et al. 2013, but see Kelly et al. 2013), which is somewhat sur¬ 
prising in view of their higher X-ray variability (e.g. Ponti et 
al. 2012). This unusually strong optical variability could in¬ 
dicate that some of the X-rays are reprocessed within a large 
scale wind. In which case the models above in Figure [8^, 
which require only disc emission to fit the optical, may be 
overestimating the mass accretion rate. Nonetheless there 
must be sufficient X-ray emission to power an additional 
optical component. The X-ray model fit based on 10®M© 
shown in Figure [8 )d has an X-ray peak which has similar 
uLv to the required far UV emission, so this model with 
L/LEdd ~ 0.56 would require that almost all the X-rays are 
reprocessed. A more likely reprocessing fraction of 10-50% 
then requires a higher Eddington ratio of 1 — 5. So perhaps 
the most plausible model is a combination of the two 10®M© 
models in Figure[7K and Figure[7)D, where the X-ray emission 
is partly suppressed by wind/advection losses and the opti- 















The Best Constraints on A Super-Eddington Accretion Flow 11 


cal emission is correspondingly enhanced by reprocessing in 
the wind. 

We note that if winds and advection do play a signifi¬ 
cant role, it is curious that the X-ray spectrum and variabil¬ 
ity properties are not more noticeably different from other 
well studied NLSl such as PG 1244-1-026, RE J1034-I-396 
and RX J0136.9-3510 (e.g. Middleton et al. 2009; Jin et al. 
2009; Jin et al. 2013). This then raises the possibility that 
these sources are likewise super-Eddington. PG 1244+026 
has an SED which is acceptably well fit by a lO^M© black 
hole accreting at about Lsdd (Done et al. 2013). This 
SED has a distinct outer disc spectral shape in the opti- 
cal/UV (e.g. Doneet al. 2013), which constrains oc 

{M^L/Lsdd)^^^■ Reducing the mass by a factor of 4 to 
2.5 X 10®M© i.e. removing the Marconi et al. (2008) ra¬ 
diation pressure correction to the H/3 line width, then re¬ 
quires L/Lsdd = 16, which clearly means that this source 
could be very super-Eddington. However, this is probably 
an over-estimate as the well defined HE break in the power 
spectrum of this source (Jin et al. 2013) suggests at least 
a factor 5 higher mass than in RX1140, but this still gives 
L/Lsdd = 4. Advection and winds could also be significant 
in this source, which would invalidate the black hole spin 
determination described in Done et al. (2013). 

If the majority of NLSl are super-Eddington, but the 
typical Broad Line Seyfert I’s have L/LEdd ~ 0.05 — 0.2, 
then this implies a distinct bimodal (or a very long tail), 
mass accretion rate distribution. Also, such large energy 
losses due to advection/winds are probably not seen in the 
majority of Ultra-luminous X-ray sources (e.g. Sutton et al. 
2015), although they may be present in the more extreme 
sources discussed by Gladstone, Done & Roberts (2009). 
The current data suggest that most of these sources rep¬ 
resent moderately super-Eddington flows onto stellar mass 
black holes with masses around SOM© compared with those 
with ~ lOM© seen in our Galaxy, due to their lower metal- 
licity resulting in less mass loss from winds: Zampieri & 
Roberts (2009). Theoretical models predict quite high ef¬ 
ficiency for these flows (Jiang, Stone & Davis 2014, but 
see also Sadowski et al. 2015 who report lower efficiency). 
Glearly there remain are many aspects of super-Eddington 
accretion flows which are not yet understood. 


7 SUMMARY 

In this paper we use new XMM-Newton data from one of 
the lowest mass AGN known to try to understand the na¬ 
ture of the soft X-ray excess, and how this fits in the context 
of its broad band SED. The mass from H/3 (uncorrected for 
radiation pressure: Marconi et al. 2008), upper limit from 
reverberation mapping and fast X-ray variability are all con¬ 
sistent with a mass of 10® M©. This, together with a steep 
2-10 keV X-ray spectrum, and even steeper soft X-ray emis¬ 
sion in this object make this a clear member of the ‘simple’ 
NLSl class of AGN, where all objects with low black hole 
mass accreting at a high accretion rate. 

We fit the time averaged X-ray data together with 
the covariance spectra of the fast variability to show that 
both low temperature Gomptonisation and highly smeared, 
ionised reflection models can fit the data, but that both re¬ 
quire an additional component at the lowest energies which 



Figure 9. relation formed by the RM sample, where 

u'^Yns Ibe excess variance of 2-10 keV for 40 ks timescale. 
values are taken from Ponti et al. (2012), the RM masses are 
taken from Peterson et al. (2004, 2005), Denney et al. (2006), Du 
et al. (2014). The black solid line is the regression line assuming 
as an independent variable. Dashed and dotted lines show 
the ±1(T and ±2 (t confidential ranges for new observation. Red 
vertical line indicates the mean for RX1140, with ilcr in¬ 

tervals determined from Obs-1 (orange) and Obs-2 (green). The 
horizontal red line indicates the mass of RX1140 predicted by the 
regression line. 

we interpret as being from the inner disc. However, the RMS, 
lag and coherence spectra all point to a break in variability 
properties between the soft and hard X-ray bands, which 
is not consistent with the reflection models. Instead these 
support the low temperature Gomptonisation model as this 
does not extend into the 2-10 keV band. All these proper¬ 
ties are similar to other ‘simple’ NLSl such as PG 1244+026 
(Jin et al. 2013), RE J1934+396 (Middleton et al. 2009) and 
RX J0136.9-3510 (Jin et al. 2009). 

However, putting the X-ray data into the context of 
the broadband SED reveals something quite new. The far 
UV and variable optical flux require an outer disc which 
has L/LEdd ~ 10 for an assumed mass of 10® M©. This far 
over-predicts the soft X-ray emission for standard efficiency 
even assuming zero spin for the black hole. This strongly 
requires substantial energy loss via advection and/or winds 
unless the black hole mass is underestimated by a factor of 
10, or the variable optical flux is not predominantly from the 
outer accretion disc (both of which seen extremely unlikely). 
If so, the other similar NLSl are also likely to be similarly 
super-Eddington. This removes the requirement for low spin 
in these objects (Done et al. 2013), but also makes it highly 
unlikely that the geometry of these NLSl are well repre¬ 
sented by a flat disc, as used in reverberation studies. 
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